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Abstract 

We deal with the periodic boundary value problem for a second-order nonlinear 
ODE which includes the case of the Nagumo type equation v xx — gv+n(x)F(v) = 0, 
previously considered by Grindrod and Sleeman and by Chen and Bell in the study 
of the model of a nerve fiber with excitable spines. In a recent work we proved 
a result of nonexistence of nontrivial solutions as well as a result of existence of 
two positive solutions, the different situations depending by a threshold parameter 
related to the integral of the weight function n(x). Here we show that the number 
of positive periodic solutions may be very large for some special choices of a (large) 
weight n. We also obtain the existence of subharmonic solutions of any order. The 
proofs are based on the Poincare - Bikhoff fixed point theorem. 
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1 Introduction 



This paper deals with the study of periodic solutions to a class of second order 
ordinary differential equations arising in the search of stationary solutions of 
a partial differential system modelling a nerve fiber with excitable spines (see 
[3], [10]). 



Preprint 



The model equation considered in [3, pp. 391-395] takes the form 

Vxx - gv + n(x)F(v) = (1) 

where g > is a given constant, n(-) is a positive /5-periodic piecewise con- 
stant function and F : R — > K. is a smooth mapping with — F having a 
iV-shape. Typically, F(s) has three zeros < a < 1 and F(s) > for 
s e (-oo,0[U]a, l[,F(s) < for s e ]0, a[U]l, +oo). According to [3, pp.391- 
392], a possible example for F is given by F(s) = f(R(s)), where f(u) = 
u(l — u)(u — a) and R is a smooth monotone increasing function with R(0) = 0, 
/2(a) = a and i?(l) = 1. Actually, equation (1) comes from the system 

v xx - gv + n(x)f(v) = 

< 

^ u = S(u;p) :=u- 

by assuming the function S'(u) invertible with inverse R(v). Such an assump- 
tion is satisfied when the spine stem resistance Rs = p' 1 approaches zero (cf. 
[3, pp.386-387] ). 

Equations of the form of (1) present a typical threshold phenomenon that may 
be easily described if we think for a moment at the weight n(x) as a constant 
function and look for constant solutions. It is clear that if n = is a small 
positive number, then the only /9-periodic solution is the trivial one. On the 
other hand, if n — rt\ is sufficiently large, the line y = gs intersects the V-- 
shaped curve y = niF(s) in two nontrivial points (at least) and therefore we 
have the existence of at least two positive /9-periodic solutions (see Figure 1). 



In [3] the authors, using a phase-plane analysis and gluing together the solu- 
tions found for different values of the coefficients, showed that a similar result 
can be obtained for a piecewise constant weight n(x) which, in a period, takes 
two values rt\ > n > (with n(x) = n on ]a, (5\ and n(x) = rt\ on ]0, a[ ). 
Then, they obtained [3, Lemma 4.1] the existence of the only trivial solution 
if a and no are small (which means that the weight is close to a small value no 
for the main part of the time) and the existence of one positive /5-periodic so- 
lution if rii is sufficiently large and a is not too close to zero (which means that 
the weight is large for an adequate amount of time). In a further result of the 
same paper, the authors [3, Lemma 4.3] also claimed the existence of "many" 
periodic solutions when S(u) is not invertible and n\ is large enough. The 
proof, however, in this case is only suggested "by reading the superimposed 
phase portraits" . For all these results, some further technical conditions (that 
we do not recall here) were assumed. Among such conditions, an hypothesis 
which implies 

/ F(s)ds>0, (2) 
Jo 



2 



Fig. 1. Intersections of y = nF(s) with y = gs. 



was required. 

A few questions may arise from these conclusions in [3], in particular, con- 
cerning the possibility of extending Chen and Bell results to a general class 
of weight functions and to the existence of at least two positive /9-periodic 
solutions when the weight is sufficiently large (in some sense), as well as to 
prove in a more general way (i.e., using some argument which is not based 
by reading the superimposed phase portraits and therefore applicable only 
to weights which are two-step functions or small perturbations of two-steps 
functions [3, Remark 4.6]) the existence of "many" periodic solutions for large 
but arbitrarily shaped profiles. 

In a previous recent article [16], we addressed our attention to the first two 
questions, as well as to some related problems. In particular, we showed that 
the same features are preserved for a broad class of functions F and for an 
arbitrary positive /9-periodic weight n G L 1 ([0, / 3]). In this case, the L 1 -norm 
(or the mean value) of n(x) plays the role of a threshold parameter in the 
sense that we have only the trivial solution when \n\i is small and at least 
two positive /5-periodic solutions when \n\i is sufficiently large. Indeed, these 
results are proved in [16] for rather general second order equations including 
(1) as a particular case. The mathematical tools employed in [16] were based 
on upper and lower solutions, the study of quadratic forms in Hilbert spaces 
and critical point theory. Also in [16] condition (2) was needed for the proof 
of the existence of two nontrivial solutions. 
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It is the aim of this work to consider now the third question, that is, to prove 
the existence of a large number of periodic solutions (both harmonics and sub- 
harmonics), for a general (i.e., not necessarily piecewise constant) positive 
periodic weight function. Note that now the elementary observation made 
above which suggested to look at the intersections of the graph of y — nF(s) 
with the line y = gs in order to prove the existence of at least two periodic 
solutions also for a non-constant weight does not work anymore. Here we have 
better to look at the phase-plane portrait of the planar system 

v ' = y 

(3) 

y' = gv — n{x)F[v) 

and observe that if n(x) = n is a constant function, with n sufficiently large and 
F\a) > then, near the point (a, 0) an equilibrium point (a„, 0) appears and 
such an equilibrium point is a local center surrounded by a family of periodic 
orbits contained in the strip ]0, 1[ xf. The fundamental period of these orbits 
becomes larger as their energy increases, a situation which is reminiscent to the 
one encountered in the study of the trajectories approaching the separatrices 
in the nonlinear simple pendulum equation (see Figure 2). 



This point of view addresses our investigation toward the search of periodic 
trajectories "near" (a, 0) also when n(x) is not constant, provided that n(x) 
has a large principal part. Indeed, we will show that such periodic trajecto- 
ries actually exist (see Theorem 3.1 of Section 3). More precisely, using the 
Poincare - Birkhoff fixed point theorem we shall prove that, if we split 

n(x) = n + n(x), (4) 

where n a suitably chosen constant value (for instance, in some situations, 
1 ff 3 

we could take n — — / n{x) dx, but other possible choices may be suitable 
p Jo 

as well, depending on the weight function) then, the number of positive (3- 
periodic solutions becomes large if n grows to infinity and |n|i is small (in 
a suitable sense). Moreover, the same fact is true also with respect to the 
subharmonic solutions (that is the m/9-periodic solutions having mf3 as their 
minimal period). To this aim, we shall consider Eq. (1) as a modification of 
the autonomous equation 

v X x - gv + nF(v) = 

and treat the remaining term n(x)F(v) as a perturbation (see, for instance, 
[2], [8], [9], [11] for a similar approach in the study of some different equations). 
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Fig. 2. Phase-portrait of system (3) for n(x) = n = 20 (constant), g = 0.1 and 
F(x) = x(l — x)(x — a), with a = 0.6. Plotting the orbit-segments at different 
initial points (xo,Q) with xq £]0, 1], but along a same time interval [— T, T], one can 
see the presence of periodic trajectories around the stable equilibrium point (a n ,0) 
which is near to (a, 0). The period needed to complete one turn becomes larger as 
the trajectories approach the separatrix. 

As a final remark, we observe that the only crucial assumption for our main 
multiplicity result of Section 3 is the existence of a point a £ ]0, 1 [ where 

F(a) = 0, F\a) > 0, 

an hypothesis which is always satisfied for the typical iV-functions considered 
in the literature. On the other hand, with our approach, we do not need any 
condition like (2) on Jq F(s) ds. As pointed out before, we recall that the 
positivity of the integral in (2) was required both in the phase-plane analysis 
approach of Chen and Bell [3] (based on the superposition of the phase por- 
traits for n = n and n = rii) as well as in the variational approach of our 
recent work (based on Ambrosetti - Rabinowitz mountain pass theorem [1]). 

Hopefully, the study of equation (1) or its variants may be relevant for the 
attempt of better understanding the transmission of the impulses along the 
nerve fibers also with respect to the fact that there are some serious pathologies 
which may arise as a consequence of a deficient distribution of myelin along 
the nerve axon (disorders of myelination). The weight function n(x) in the 
Nagumo equation represents the axial distribution of the myelin along the 
axoplasm. It is known that some "fat" areas alternate with some gaps (the 
so-called nodes of Ranvier). Such a profile for the distribution of n(x) has lead 
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some authors (like in [3], [10]) to the study of a piecewise constant weight 
which alternates between a small and a large value. However, the real images 
of a myelinated axon show that the true shape of the profile sometimes may 
be quite far from this idealized picture and therefore results like those in [3] 
could not be applied. Our theorems show that those conclusions are still valid 
for general positive weights and Wshaped functions. 

2 Preliminary results and notation 

We recall in this section an auxiliary result (taken from [16]), about the peri- 
odic solutions of (1), which is useful for what follows. 

The nerve fiber equation (1), has the form of a second-order equation 

v" + h(x, v) = (5) 

where /i : f xl ^lisa Caratheodory function which is /9-periodic in its first 
variable, that is, h(-, s) is measurable for all s G R and there is (3 > such that 
h(x + (3, s) = h(x, s) for almost every i6R and all s£l, h(x, •) is continuous 
for almost every x G [0, f3] and, for every r > there is a measurable function 
p r G L 1 ([0, R + ) such that \h(x,s)\ < p r (x) for almost every x G [0, f3] and 
every s G [— r, r\. Solutions of (5) are considered in the Caratheodory sense 
too. 

We look for the existence of m/?-periodic solutions to Eq. (5), for some positive 
integer m, that is for solutions of problem 

(P m ) v" + h(x,v) = 0, v(x + m(3) = v(x), Vrc G R, 

or, equivalently, v" + h(x,v) = 0, with v(m(3) — v(0) = v'(mf3) — v'(0) = 0. 

Clearly, in the model of our interest, we have 

h(t,s) := -gs + n(t)F(s) (6) 

and the iV-shape of —F(s) implies that for h(t, s) the following conditions are 
satisfied. 

(Ai) h(x,0) = 

and, for a.e. x G [0, 0\, 

(A 2 ) h(x,s) > 0, Vs < 0, h(x,s) < 0, Vs > 1. 

Under these assumptions, we obtain: 
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Lemma 2.1 [16] Suppose (Ai) and(A2) hold. Then, any possible solution v(-) 
of problem (P m ) satisfies < v(x) < 1, Vx G M.. Moreover, if h is (locally) 
lipschitzean at s = 0, then any nontrivial solution of (P m ) is strictly positive 
and, if h is (locally) lipschitzean at s = 1, then any solution of (P m ) is strictly 
less than one. 

The proof is based on some direct estimates and classical arguments from the 
theory of upper and lower solutions (see, e.g., [4], [5]). We refer to [16] for 
more details and remarks. 

Thanks to Lemma 2.1 we can confine ourselves to the solutions v(x) belonging 
to the interval [0, 1] and, moreover, we can modify as we like the function 
h(t,s) (respectively the function F(s) in (1) ) on s G M. \ [0,1] and such a 
modification will have no effect on the existence of the periodic solutions as 
long as the sign conditions outside the interval [0, 1] are preserved. 

In the sequel, standard notation is used. We only warn that for a /9-periodic 
locally integrable function «(•), we denote by its L 1 -norm on a interval of 
length (3. Sometimes (when no confusion may occur) the same symbol will be 
used to denote the L 1 -norm on [0, m/3], when we are interested in the search 
of m/5-periodic functions. 



3 Multiplicity results 



Now we focus our attention on equation (1) and assume that g > is a fixed 
constant and F : R — > H. is a continuously differentiable function such that 
F(0) = F(a) = F(l) = with a : < a < 1 and such that F(s) > for 
s E (— oo, 0[U]a, 1[ , F(s) < for s G ]0, a[ U ]1, +oo). We also suppose that 
n(-) : M. — > M + is a /3-periodic locally integrable function. 

We are going to show an application of the Poincare - Birkhoff fixed point 
theorem to the search of the periodic solutions of (1). To this aim and for 
technical reasons it will be convenient to modify suitably the function F out- 
side the interval [0, 1]. Precisely, let us set 

S(s) := max{0, min{s, 1}} 

and define 

F (s) :=F(5(s)) + k £(s), 
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where £(s) is given by 



£(s) = 



exp(l/s), s < 

0, < s < 1 (7) 

-exp(l/(l -s)), s > 1, 



and 



< kn < c := max \F(s) 
se[o,i] 



With these positions, we have that F : R — > R is a bounded (globally) 
Lipschitz continuous function which is smooth in ]0, 1[, Fq(s) = F(s) for all 
s E [0, 1], F (s) > for s < and F (s) < for s > 1. Moreover, 

sup |F„(s)| = max |F (s)| = max |F(s)| = c . 
seM se[o,i] 

A straightforward application of Lemma 2.1 ensures that for each positive 
integer m, all the possible nontrivial m/5-periodic solutions of equation 

v" - gv + n(x)F (v) = (8) 

have range into the open interval ]0, 1[ and so they are, indeed, solutions of 
equation (1). 

The fact that after our modification the nonlinearity is now bounded, suggests 
the possibility of splitting n(x) as in (4) and writing equation (8) as 

v" — gv + nF (v) = p(x, v) (9) 

which looks like a perturbation of the autonomous equation 

v" - gv + nF (v) = (10) 

by an external source 

p(x,v) := — n(x)F (v). 
Now we are in position to prove our main multiplicity result. 

Theorem 3.1 Assume that 

F'(a) > 0. 

Then, for every integer N > 1, there exists a (large) value fj, — (j,n > such 
that for each fj,2 > A*i > A* there is a (small) value e = £at iM1iM2 > such that 
for each n(-) with 

n G [yUi,/i2] and \h\\ < e 

there are at least 2N solutions of equation (1) which are j3 -periodic and take 
values in ]0, 1[ . 
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The proof of Theorem 3.1 is based on the Poincare - Birkhoff fixed point 
theorem and will be split into some steps. 



First of all, we write equation (8) as a first order planar system of the form 

x — y 

(11) 

% y = -h(t,x) 

with 

dx 

h(t,s) = -gs + n(t)F (s), x = — . 

Note that we have changed the name to the variables, by denoting the x and 
the v variables in (8) as t and x, respectively. According to our hypotheses, 
/i:MxIR^Misa Caratheodory function which is /9-periodic in t and satisfies 
a global Lipschitz condition 

\h(t, si) — h(t, s 2 )\ < A(t)\s 1 — s 2 |, for all s 1 ,s 2 G R and for a.e. t G R, 

where A(-) is a suitable /9-periodic measurable function with A(-) G L 1 ([0, / 9]). 
In our case, we actually have A(t) = g+n(t)L where L is a Lipschitz constant 
for F (which is globally lipschitzean on R). As a consequence, we have that for 
each zq = (x ,y ) G M 2 and t G R, system (11) has a unique solution ((t) = 
C(t,to,zo) = (x(t, to, zo), y(t, to, zq)) satisfying the initial condition C(*o) — z o , 
with ((t) defined for all t G R. Hence, the Poincare's map 

: Zq (->• C(/3, 0, Z ) 

is well defined as a homeomorphism of IR 2 onto itself. 

Next, we introduce a rotation number associated to with respect to a given 
point q G M 2 . 

Let m G N be a positive integer and let go = (Qi, Q2) an( i z o — ( x o, Ho) be such 
that 

C(t,0,z o ) ^q Vte[0,mP\. 
We define the rotation number (see, for instance, [7], [15]) as 

1 p 3 ? /(t) 2 + 3 ;(t)Mt, 3 ;(t)) 
rot m (, ,go)-^/ o m __ g o )2 - m _ W dt 

which represents the normalized angular displacement around the point q of 
the solution ((t, 0, Zq) for t varying along the time interval [0, m/3]. In fact, if we 
use the Priifer transformation and write the solution ((t) in polar coordinates 
(9, p) with respect to the point q , we have that p{t) 2 = \ \((t,0,z ) — q \\ 2 = 
(x(t, 0, z ) - g?) 2 + (y(t, 0, z ) - q^f > for all t G [0, m/3] and therefore the 
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number 9(t) — 9(0) is well defined. It turns out that 

rot m (z ,q ) = — ®( m ^ 
and this number counts as positive the clockwise rotations around the point 

go- 

Our main tool to prove the existence of periodic solutions is the following 
result which is adapted to our situation from W. Ding's generalized version of 
the Poincare - Birkhoff theorem [8]. 

Theorem 3.2 Let q Q = q%) G M 2 and let D be an open neighborhood of 
q such that 

C(Mo,*o)^5o, VzoEdDo, Vt o e[0,mP[,Vte[t o ,mP]. (12) 

Suppose that T C M. 2 \ D is a simple closed curve which is star-shaped with 
respect to q Q and there exist a positive integer j such that 

rot m (zi,g ) > j, Vzi G T. (13) 

Furthermore, let us assume there is a (sufficiently large) radius R > 0, with 
r C B(q ,R) such that 

Tot m (z 2 ,q ) < 1, V 1 1^2 1 1 = R. (14) 

Define A to be the open annulus bounded by V and dB(q , R). Then, there are 
at least 2j solutions (k and (t (for k — 1, . . . ,j) of (11) which are mf3-periodic 
and such that Cfc(0),Cfe(0) e A and 

rot m (2 ,g ) = k, for z = ( k (0) or z = 4(0). 

In the particular case when q = (??,0) we have also that, setting Cfc(0 = 
(uk(t),yk(t)) and (kit) = (v,k(t),yk(t)), it follows that u k and u k are mf3- 
periodic solutions of equation 

u" + h(t, u) = 

with itkit) — <2i and u k (t) — q\ having precisely 2k simple zeros in the interval 
[0,mP[. 

Proof. Let : 1R 2 — > R 2 be the Poincare's map associated to system (11). 
By the Liouville theorem it follows that <fi is area-preserving and therefore 
<p m is an area-preserving planar homeomorphism as well. Let E be the open 
bounded set confined by the curve T. By the assumption, qo G Do C E 
and OEq = T. Note also that the rotation number rot m (,2o,go) ^ s we ^ defined 
for each z G M 2 \ E due to the fact ((t,0,z ) ^ q for all t G [0, mf3] 
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and each zq $l Eq . Indeed, if, by contradiction, there is zq ^ Eq such that 
£(t*) = ((t*,0,z ) = q for some t* e]0,m/3], then C(ti,0,z ) = *i G d-D 
for some ^ G [0, £*[. This means that C(^*) = C(**j *ij z i) — Qo and we have 
a contradiction to hypothesis (12). Now, if we restrict ip := <p m to the closed 
disk B[q ,R] we have that 

and ^ satisfies the twist condition with respect to the inner and the outer 
boundaries of the closed annulus A for each integer k e [1, j}. Namely, we 
have rot m (-,g ) > k, on V and rot m (-,g ) < k, on dB(q ,R). The W. Ding's 
generalized version of the twist theorem guarantees the existence of at least 
two fixed points in A for if). These fixed points are initial values (at the time 
t = 0) of two m/9-periodic solutions Cfc(') and Cfe(') °f system (11), respectively. 
Moreover, the rotation number associated to (k an d Cfc is equal to k. 

In the special case when q = (?°, 0), and by virtue of the particular form of 
system (11), we know that the result about the rotation numbers implies that 
the first coordinates Uk and Uk of (k and (k crosses the line u = q® exactly 2k 
times in the interval [0, mf3 \ . □ 

Remark 3.1 Clearly, if u(-) is a m/3-periodic solution of the equation u" + 
h(t, u) = with h(t + f3 : u) — h(t, u) then 

u(- + 0), . . . ,u{- + jP), . . . ,u(- + (m - 

are m/3-periodic solutions as well. We consider these solutions as equivalent 
each other and we say that they belong to the same periodicity class. A further 
consequence of the Poincare - Birkhoff theorem (as pointed out in [14]) ensures 
that the solutions Uk and iik (for which the existence is claimed in Theorem 
3.2 do not belong to the same periodicity class. Obviously, also the solutions 
Uk , Uk ue , ue for £ ^ k belong to different periodicity classes (in fact their 
rotation numbers are different). We refer also to [13], [18] for a throughout 
discussion concerning this aspect. 

Moreover, we observe that the information about the associated rotation num- 
bers permits obtain some conclusions about the minimality of the period. For 
instance, if m > 2 and k > 1 are co-prime numbers (that is m/k is not further 
reducible), then it is possible to prove that the m/5-periodic solutions Uk and 
Uk are not i/5-periodic for each i — 1, . . . ,m — 1. In particular, for k = 1, we 
have that the solutions we find have m/3 as their minimal period (see [7] where 
this discussion is carried on with more details). 

In the proof of Theorem 3.1 using Theorem 3.2 we take 

g = (a fi , 0) 
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with < a n < 1 such that 

-Pb(fln) _ 9_ 
a n n 

Note that 



a n > a, Vn > and lim a n = a. 



Our next result shows that large solutions rotate slowly. 

Lemma 3.1 There is R = Ro(m, \n\i) > such that for each initial point z 2 
with 1 \z 2 \ 1 > Ro 

rot m (z 2 ,q ) < 1, V||g ||<l- 

Proof. First of all, we recall a well known consequence of the global existence 
of the solutions (cf. [12]), that is, there is a continuous nondecreasing function 
r] : IR + — > IRq , with r)(r) > r for all r > 0, snc/i £/ia£ 

||C(Mo,3b)|| > i?, Vt,t e [0,m/3], V jj^oll > ^(^)- 

The function 77 depends upon the constants which bound the growth of h(t, x) 
and hence, we have r\ = ?7 m ,|n|i (see also [17, Lemma 2]). 

Let R > be such that ||C(t)|| > R for alU e [0,m/3]. We claim that if R > 1 
then 

rot m (C(0),0) < 1/2. 

Indeed, let us assume the contrary and suppose that rot m (C(0), 0) > 1/2. This 
implies that the projection C(0/llC(0ll °f C(0 to S* 1 covers at least one of the 
intersections of S 1 with the standard quadrants in the xy-pl&ne. Just to fix 
the ideas, let us assume that the trajectory crosses the first quadrant. Hence, 
there are < t\ < t 2 < m(5 such that 

C(ti) = {xih)^)) = (l, yi ), with yi > (R 2 - 1)5 
C(t 2 ) = (x(t 2 ),y(t 2 )) = (x 2 , 0), with x 2 > R 
\\C(t)\\>R, x(t)>l, y(t)>0, Wte]h,t 2 [. 

Then, passing to the polar coordinates (via the Priifer transformation), we 
find, for t E [ii, t 2 ] and using the fact that h(t, s)s < 0, for all s > 1, 

_ y(t) 2 + x(t)h(t,x(t)) yjtf _ 

{ ) ~ xW+vW - x(ty + y (ty ~ 8111 

For t £ [ti,t 2 [, we have 
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cotflt )-cot(0(fi )=/ —s- = - / -^frr ds 

Je(t) sm 2 u Jti sin^(s)) 



< f lds = t-t 1 < m(3. 
Jti 

Therefore, we get 

cot(0(f)) - — 1 -r < cot(0(*)) - yr 1 < Vt G M 2 [. 

(it" 1 — 1)2 

This yields to a contradiction as t — > and therefore our claim is proved. 

Adapting a result in [7, Lemma 2.2] to our situation, we have that for every 
e > there is _R £ > 1 such that 

|rot m (C(0),0)-rot m (C(0),g )| <e, 

holds for every point q = (%,0) with ||g || < 1 and all solutions ((■) with 

\\C(t)\\>R e ,Vte[o,mf3}. 

Thus there exists a sufficiently large radius i2* > 1 such that 

rot m (C(0),g ) < \, V||g || < l,VCwith ||C(t)|| > R* , Vt G [0,m/3] 

and, by the result recalled at the beginning of the proof, we achieve the con- 
clusion by taking 

Ro > r](R*). 

Thus the proof is complete. □ 

As a next step, we evaluate the rotation numbers of small solutions around 
Qo = ( a n, 0). To this end, we consider autonomous equation (10) that we write 
as a first order system 

x — v 

(15) 

y = gx- nF (x). 
System (15) is a conservative system with energy 

11 f x 

E n (x, y) = -y 2 - -gx 2 + fUFo(x), with F (x) = / F (s) ds. 
II Jo 

For simplicity in the notation, in the sequel we set E n = E. Let us fix d > 
and b G ]a, 1[ such that 

Fq(x) > d , Vi G [a, b] 
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and take 

. g 

do 

such that 

a < a n < b, for n > /x • 

For n > yU we have that £"(-,0) is strictly increasing on [a, b] and vanishes 
in a n . Hence, on the interval [a,b] we have that E(-,0) is a strictly convex 
function having absolute minimum at the point a n . If we take now a constant 
c n with 

E(a n ,0) <c n < mm{E(a, 0), E(b, 0)}, 
we have that the level line 

T n = {(x,y) : a<x < b, E(x,y) = c n } 

is a star-shaped curve around the equilibrium point go • Moreover, for any 
point g G T n we have that the solution of system (15) with (x(0),y(0)) = g is 
periodic and the corresponding orbit coincides with T n . 

In order to avoid misunderstanding with the previously defined rotation num- 
ber, we denote by Mot m the rotation number associated to (10) along the time 
interval [0,m{3]. If we denote by r n the fundamental (i.e., minimal) period of 
Y n , then we can conclude that 



&ot m (q,q ) > 



mf3 



Vg g r. 



We claim that lim r n = 0. 

Indeed, by the well known time-mapping formula we have 

du r- f b t du 



where 



V2 r — * — r+ vi i 

JbZ (c*-FAu. 0^2 Ja 



^ (Cs-f?(«,0))5 ^ (Cs-£?(U,0))3 



&r 6+ E(b^,0) = E(bt,0)=c n . 

For -u G [a s , we have 

c s - f?(w, 0) = £?(&+, 0) - E(u, 0) = / E'(x, 0) da; 

= f bn {E'{x,0) -E'{a n ,0))dx = /"*" ( [* E"(£,0)dA dx 
> ^f 1 ((# " an? ~ (u - a,?) . 
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Therefore, via an elementary integration, we obtain 



7T 



tT< 

Vnd - g 



A similar computation for r n yields to 



2tt 

n ~ y/nd - g 



and this proves the claim. Hence we have that 

mf3\JndQ — g 



Mot m (q,q ) > 



2tt 



Wq eT n , 



which shows that £%ot m (q, go) — * +00 with the order of \fn. 



Now we are in position to complete the proof of Theorem 3.1. 

Fix m > 1 and N > 1. Let n = hn > /-to be such that 

@ot m (q, g ) = crn > a > N, Vg G T n , 

holds for each n > jj,. Fix also an interval 112] c]/i, +00) and consider 
equation (9) with n G [//i,/^]- By the continuous dependence of the solutions 
from the coefficients and from initial data (otherwise, some direct estimates 
may also be performed) we have that there exists E\ > (depending on jii and 
/i 2 and, in turns, also on N) such that for each forcing term p with \p\i < E\ , 
condition (12) is satisfied as well as rot m (zo, go) > N holds for each g = (a^, 0) 
and zo G T n provided that n G [/ii, 112] - Recalling the definition of Co as a bound 
for |F |, we have that for 

i~i / £l 
\n 1 < e = — , 

it follows that (12) holds and 

rot m (zi, g ) > N, V^i G T n . 

For any chosen n(x) with n G [^1,^2] and |n|i < 5, we can take a sufficiently 
large radius R n such that T n C S(g , i? n ) and 

rot m (z 2 ,g ) < 1, V||z 2 || = -Rn- 

The Poincare - Birkhoff theorem (Theorem 3.2) ensures the existence of 2N 
solutions which are m/3-periodic. Clearly, for m — 1 we have exactly the result 
claimed in Theorem 3.1 and thus the proof is complete. □ 

Clearly, from the above proof we have a result about the existence of sub- 
harmonic solutions (of period m(5 with (5 > 2) as well. Indeed, the discussion 
about the minimality of the period given in Remark 3.1 yields to: 
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Theorem 3.3 Assume that 

F'{a) > 0. 

Then, for every integer m > 2 and each K > 1, there is a (large) value fi = 
Hk > such that for each /i 2 > A*i > A* ^ere a (small) value e = Sk,^,^ > 
snc/i t/iat /or each n(-) 

n G fj, 2 ] and \h\i < e 

there are at least 2K solutions of equation (1) which are mf3 -periodic and take 
values in ]0, 1[ . Moreover, such solutions do not belong to the same periodicity 
class and also they are not if3-periodic, for each i = 1,2, ... ,m — 1, so that 
the period mfi is minimal. 

Proof. Once we have fixed m and K, we take the set 

N ;= {I, = 1< l 2 < . • . < l K ;= N} 

made by the first K numbers which are co-prime with m. 1 According to 
the conclusion of the proof of Theorem 3.1, we find that there are 2N > 2K 
periodic solutions of period mf3. As a consequence of Remark 3.1 we know that 
these 2N solutions appears in pairs which are characterized by the fact that 
they have, respectively, 2, 4, 6, ... , 2N oscillations in the time interval [0, m(5[ . 
Such oscillations are associated to the rotation numbers of the solutions (an 
information which comes from the use of the Poincare - Birkhoff fixed point 
theorem). As we already explained in Remark 3.1 we are sure that for those 
solutions having a rotation number I which is co-prime with m, it holds that 
m(3 is the minimal period. By the choice of iV and the set M we are lead 
to conclude that at least 2K among the 2N periodic solutions (that is those 
with rotation numbers equal to l x , l 2 , . . . , Ik) are of minimal period. This 
observation concludes our proof. □ 

We end our paper, by showing that it is possible to realize a "good" splitting 
of the weight function n(x) for a profile like the one considered in [3]. 

Example 3.1 Consider now equation (1) and assume, like in [3] that n(x) is 
a periodic piecewise constant function satisfying 

n(x) = 

with < no < n\ . Like in [3] we assume that F is a smooth function satis- 
fying the sign conditions we have already described at the beginning, that is 
F(0) = F(a) = F(l) = with a : < a < 1 and such that F(s) > for 

1 Two positive integers I and m are co-prime (or relatively prime) if GCD(/, m) = 1 



n\ if x G ]0, a[ mod (5 
no if x G ]a, (5\ mod (3 



(16) 
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s G (— oo, 0[ U ]a, 1[ , F(s) < for s G ]0, a\ U ]1, +00). In order to apply our 
result, we also suppose that F'(a) > (an hypothesis which is always satis- 
fied by the nerve fiber models considered in [3], [10]). In [16] we have already 
proved the existence of threshold values, associated to a general shape of n(x) 
which imply, for the particular case of (16), that the only periodic solution is 
the trivial one when no and a are sufficiently small, while we proved the ex- 
istence of at least two nontrivial /5-periodic solutions (under a weak technical 
assumption on F(s)) when ri\ is sufficiently large and a is close to (3. Now, in 
view of Theorem 3.1 and Theorem 3.3 we propose the splitting 

n(x) — ni + n(x) 

where 

!0 if x G]0,a[ mod (3 
n — rii iixE\a,(3[ mod j3 
By this choice n = ri\ and, moreover, we have 

Ni = \nW([o, m p\) = m ( n i - n o)(f3 - ot). 

In spite of the fact that we take n — n\ large, we are allowed to make \h\\ as 
small as we like, by taking a sufficiently close to f3. Hence Theorem 3.1 and 
Theorem 3.3 can be applied and the existence of a large number of harmonic 
and "true" subharmonic solutions for the equation (1) is guaranteed. 
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